On the origin of the non-monotonic doping dependence of the in-plane resistivity 
anisotropy of Ba(Fei_j.T^)2As2, T = Co, Ni and Cu 
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The in-plane resistivity anisotropy has been measured for detwinned single crystals of 
Ba(Fei-a;Ni2,)2 As2 and Ba(Fei_iCua;)2As2. The data reveal a non-monotonic doping dependence, 
similar to previous observations for Ba(Fei_a;Coa;)2As2. Magnetotransport measurements of the 
parent compound reveal a non-linear Hall coefficient and a large linear term in the transverse 
magnetoresistance. Both effects are rapidly suppressed with chemical substitution over a similar 
compositional range as the onset of the large in-plane resistivity anisotropy. This suggests that the 
relatively small in-plane anisotropy of the parent compound in the spin density wave state is due to 
the presence of an isotropic, high mobility pocket of the reconstructed Fermi surface. Progressive 
suppression of the contribution to the conductivity arising from this isotropic pocket with chemical 
substitution eventually reveals the underlying in-plane anisotropy associated with the remaining 
Fermi surface pockets. 

PACS numbers; 74.25.F-, 74.25.fc, 74.70.Xa, 75.47.-m 



o 

> 

in 
m 

rn 
O 



X 



I. INTRODUCTION 

Recent measurements of detwinned single crystals of 
Ba(Fei_a;Coa;)2As2 revealed a non-monotonic doping de- 
pendence of the in-plane resistivity anisotropy In 
striking contrast, the lattice orthorhombicity diminishes 
monotonically with increasing Co concentration^, rais- 
ing the question of the origin of the nonmonotonic be- 
havior of the in-plane resistivity anisotropy, and the ex- 
tent to which it is, or is not, generic to this family of 
compounds. Here we present measurements of the in- 
plane resistivity anisotropy of the closely related cases of 
Ni and Cu-substituted BaFe2As2 which reveal a similar 
non-monotonic compositional dependence of the resistiv- 
ity anisotropy. Furthermore, magnetotransport measure- 
ments indicate that the effect is closely coupled to the 
progressive erosion of the contribution to the conductiv- 
ity from an isotropic, high- mobility, reconstructed Fermi 
surface (FS) pocket. 

BaFe2As2 is a representative "parent" phase of the Fe- 
pnictide superconductors^^—. The material has an anti- 
ferromagnetic (AFM) ground state, comprising stripes 
of ferromagnetically aligned moments which alternate 
antiferromagnetically along the orthorhombic a-axis^. 
Back folding of the bands according to the antiferromag- 
netic wavevector results in a reconstructed Fermi sur- 
face consisting of several small pockets, as evidenced 
by both Angle Resolved PhotoEmission Spectroscopy 
(ARPES)^ and quantum oscillations^. The material 
has a Neel temperature T/v close to 140 K, with val- 
ues depending slightly on growth conditions and an- 



nealing treatments^iiSr— . Significantly, the Neel tran- 
sition in BaFe2As2 is accompanied by a tetragonal-to- 
orthorhombic structural transition®. For the specific 
cases of Co, Ni and Cu substitution relevant to the cur- 
rent work, the structural transition occurs at a slightly 
higher temperature Tg than the magnetic transition, with 
a temperature difference that monotonically increases 
with increasing concentration of the substituent, at least 
until the top of the superconducting dom o-'^^'^^d^" — . The 
origin of the splitting of Tg and T^r with chemical sub- 
stitution is not clear, but consideration of the effect 
of crystal quality on the splitting of the transitions in 
CeFeAsOi^ implies that this effect might, at least in 
part, be associated with the strong in-plane disorder in- 
troduced by partial substitution on the Fe site. In both 
families, the structural transition breaks a discrete ro- 
tational symmetry (C4 to C2) of the high-temperature 
phase without introducing a new translational symme- 
try, and is widely referred to as a nematic transition, 
borrowing language from the field of liquid crystala^^. 
Understanding the origin of this effect is a key compo- 
nent of a complete theoretical description of the occur- 
rence of superconductivity in this family of compounds, 
motivating both theoreticaP^"^^ and experimental"— 
investigation of the nematic transition and the associated 
in-plane anisotropy. 

BaFe2As2 tends to form dense structural twins on cool- 
ing through Ts, corresponding to alternation of the or- 
thorhombic a and b axes through the crystal^. The rel- 
ative twin population can be influenced by application 
of an in-plane magnetic field due to the in-plane suscep- 
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tibility anisotropy associated with the colinear antiferro- 
magnetic structure"^ ^. However, the degree of detwinning 
that can be achieved for typical laboratory fields is only 
modest and the anisotropy can only be explored for tem- 
peratures below Tpf^^. Much larger changes in the rela- 
tive twin domain population can be achieved by use of 
uniaxial mechanical stress, which also permits measure- 
ment of the resistivity anisotropy through Tj^^^. Such 
measurements reveal a relatively small in-plane resistiv- 
ity anisotropy in the parent compound BaFe2As2^ i^^'^^ , 
with the resistivity along the ferromagnetic direction pt, 
slightly greater than that along the antiferromagnetic di- 
rection Pa- In contrast, the anisotropy Pb/pa is initially 
found to increase with Co substitution, despite the fact 
that the orthorhombicity (a — b)/2{a + b) monotonically 
decreases with increasing Co concentration^'^^. As an- 
ticipated, the anisotropy eventually diminishes to unity 
when Ts is completely suppressed. Perhaps coinciden- 
tally, the maximum in pt/pa is found to occur for a Co 
concentration close to the beginning of the superconduct- 
ing "dome"—. 

Optical reflectivity measurements for mechanically de- 
twinned crystals reveal that the origin of the resistiv- 
ity anisotropy is principally caused by changes in the 
spectral weight (i.e. is due to changes in the FS mor- 
phology), rather than by changes in the scattering^^. 
The anisotropy is present over a wide energy scale, 
clearly involving whole bands rather than just the be- 
havior at the Fermi energy, and the dichroism is largest 
for the parent compound, consistent with the struc- 
tural orthorhombicity^^. In addition, recent polarized 
ARPES measurements on detwinned crystals reveal that 
the structural transition is associated with an increase 
(decrease) in the binding energy of bands with principal 
dxz (dyz) character^^'^®. The relative degree of splitting 
of the two bands is largest for the parent compound and 
decreases monotonically with increasing Co concentra- 
tion. 

Despite these recent advances, the origin and signif- 
icance of the nonmonotonic doping dependence of the 
in-plane resistivity anisotropy in Ba(Fei_a;Co2;)2As2 re- 
mains unclear. Initial calculations with a net orbital 
polarization have suggested the magnitude of resistivity 
anisotropy could strongly depend on the density of states 
near the Fermi level, which does not necessarily have a 
monotonic doping dependence^. In this regard, it is es- 
pecially useful to consider other dopants in order to es- 
tablish systematic trends. In this paper, we compare the 
temperature and doping dependence of the in-plane resis- 
tivity anisotropy of detwinned crystals of Co, Ni and Cu 
substituted BaFe2As2. All three reveal a non- monotonic 
doping dependence to pt/pa, but with the maximum in- 
plane anisotropy occurring for different ranges of the 
dopant concentration. For the two cases of Co and Ni 
substitution, for which a direct comparison of the physi- 
cal properties is best motivated both phenomenologically 
and also based on ab-initio calculations- ^ , the onset of the 
large in-plane resistivity anisotropy coincides with a sup- 
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FIG. 1: (Color online) (a) Image of the surface of an un- 
stressed single crystal of BaFe2As2 at a temperature below 
10 K, obtained using polarized-light microscopy, as described 
in the main text. Stripes of light and dark contrast are asso- 
ciated with different twin orientations. Arrows indicate the 
orientation of crystal a and 6 axes, (b) Image of an un- 
stressed crystal of Ba(Fei-2,Nij;)2As2 with a;=0.014, also re- 
vealing horizontal stripes corresponding to the two twin orien- 
tations, but with a significantly reduced contrast. The jagged 
feature running from top to bottom of the image is due to 
surface morphology. To better reveal the stripes, the inten- 
sity was integrated in the horizontal direction and plotted 
as a function of vertical position, y (right hand axis), (c) 
The same region of the same crystal as shown in panel (b), 
but with uniaxial stress applied in the direction indicated, 
revealing the detwinning effect of the uniaxial stress. Note 
the absence of horizontal stripes in both the photograph and 
also the integrated intensity plot, (d) The in-plane resistiv- 
ity anisotropy expressed as Pb/pa for three different crystals 
of Ba(Fei_a;Nia;)2As2 with x ~ 0.017, illustrating that the 
measurements are reproducible. 

pression of both a large non-linear contribution to the 
transverse resistivity pxy, and also of a linear term in 
the magnetoresistance, suggesting a common origin for 
the two effects. 



II. EXPERIMENTAL METHODS 

Single crystals of Ba(Fei_j;Nia;)2As2 and 
Ba(Fei_a:Cua:)2As2 wcrc grown from a self-flu»i2,il2i. Ba 
was combined with a mixture of FeAs and Ni/Cu with a 
ratio of Ba:FeAs:Ni/Cu — l:4:x. The mixture was held 
in an alumina crucible and sealed in quartz, and was 
heated slowly to 1190 °C and then cooled down to 1000 
°C in 60 hrs, at which temperature the remaining flux 
was decanted. The Ni and Cu content of the resulting 
crystals was measured by electron microprobe analysis 
(EMPA) using BaFe2As2 and elemental Ni and Cu as 
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standards. Measurements were made for several points 
on each crystal, with standard deviations for the Ni and 
Cu concentrations which were generally below 10% and 
14% of the absolute values respectively. 

A mechanical detwinning device, with a design sim- 
ilar to that which was previously described for our 
earlier measurements of Ba(Fei_2;Co2;)2As2i but mod- 
ified to be suitable for somewhat smaller crystals of 
Ba(Fei_:rNia;)2As2 and Ba(Fei_2;Cux)2As2, was used to 
mount the crystals for measurement of the in-plane re- 
sistivity. Crystals were cut into rectilinear bars with the 
tetragonal a-axis at 45 degrees to the sides of the bar 
(corresponding to the natural crystal facets), and with 
in-plane aspect ratios of approximately 1:1.2. The crys- 
tals were placed on a horizontal platform on the detwin- 
ning device, such that an insulated Cu plate rested on 
the edge of the crystal. Uniaxial pressure was applied by 
tightening the Cu plate against the edge of the crystal. A 
similar magnitude stress was applied to all crystals stud- 
ied, but absolute values could not be estimated for the 
small cantilevers used in this study. On cooling through 
Ts, the uniaxial pressure favors the twin orientation with 
the shorter 6-axis along the direction of the applied com- 
pressive stressi. 

The degree of detwinning was monitored for several 
representative crystals via polarized light microscopy. 
Samples were positioned on the cold finger of a vac- 
uum cryostat. The samples were illuminated with ~800 
nm light, linearly polarized at approximately 45 degrees 
to the orthorhombic a/b-axes. The reflected light was 
passed through an optical compensator and analyzed 
by an almost fully crossed polarizer to maximize the 
contrast in birefringence between the two twin orienta- 
tions. Whereas good contrast between the two twin ori- 
entations is possible for the parent compound (FiglIJa)) 
and for Ba(Fei_a;Coa;)2As2 for compositions across the 
phase diagram^l, Ni and Cu substitution are found to 
rapidly suppress the contrast between the two twin ori- 
entations. Representative images for a single crystal of 
Ba(Fei_2;Ni2:)2As2 with x — 0.014 with and without uni- 
axial stress are shown in Figjl] (b) and (c) respectively. 
The application of uniaxial stress clearly results in a sin- 
gle twin domain orientation over the field of view, in this 
case about 22 microns. The origin of the reduced con- 
trast relative to Ba(Fei_a;Coa;)2As2 is not known. The 
twin domain dimensions seen in Figs[TJb) are similar to 
those found for Ba(Fei_2:Coa;)2As2, so it is unlikely that 
the width of the domains for higher Ni concentrations 
fall below the resolution of the microscope, although we 
cannot completely exclude this possibility. Although the 
optical imaging was not possible for x > 0.014 for Ni- 
substituted samples, or for x > for Cu substituted 
samples, a clear change in the resistivity was observed 
as a consequence of the applied stress for all underdoped 
compositions. Since we are unable to independently de- 
termine the degree of detwinning for these compositions, 
the observed in-plane resistivity anisotropy must there- 
fore be considered a minimum bound on the actual value. 



However, for several representative compositions, mea- 
surements were made with multiple crystals, in each case 
revealing almost identical changes in the in-plane resis- 
tivity for stress applied parallel and perpendicular to the 
current (FiglTJd) for example), providing evidence that 
even for the samples for which optical characterization 
failing to verify the effect of detwinning, the samples were 
essentially fully detwinned. 

The in-plane resistivity was measured using a stan- 
dard 4-probe configuration. Crystals were rotated so 
that the applied pressure was parallel and perpendicu- 
lar to the current, enabling measurement of pb and pa 
respectively. The same contacts were used for each mea- 
surement, to avoid uncertainty in the geometric factors. 
Measurements were made for multiple crystals of sev- 
eral compositions to ensure reproducibility of the results. 
The in-plane resistivity was also measured for unstrained 
crystals in order to determine Tg for each composition. 
For several samples, the in-plane resistivity after strain- 
ing was also measured, and no change in either Tg or T/v 
was observed. Finally, both p^x and the transverse resis- 
tivity pxy were measured at the National High Magnetic 
Field Laboratory (NHMFL) in Tallahassee in dc mag- 
netic fields up to 35 T and in Stanford in fields up to 14 T. 
The magnetic field was always oriented along the c-axis 
and the samples were mounted using a 6-point contact 
configuration. Measurements were made for both posi- 
tive and negative field orientations in order to subtract 
any small resistive component due to contact misalign- 
ment. 



III. RESULTS 

Results of measurements of the in-plane resistivity for 
representative single crystals of Ba(Fei_j;Nia;)2As2 and 
Ba(Fei_2,Cu2;)2As2 held under an applied uniaxial stress 
in the detwinning device are shown in Figsl2ja) and (b) 
respectively. For each composition, pa and pb were mea- 
sured for the same crystal. Nevertheless, the small crys- 
tals used for these measurements are susceptible to dam- 
age while being repositioned in the detwinning device, so 
the resistivity has been normalized by its value at 300 K 
in order to avoid uncertainty arising from subtle changes 
in the geometric factors between each measurement. For 
cases for which the applied stress is perpendicular to 
the current, the data are labeled as pa (green curves), 
whereas for cases for which the applied stress is parallel 
to the current the data are label- led as pb (red curves). 
As described previously^ T/v is unaffected by the small 
stress used to detwin the crystals, but the structural tran- 
sition is rapidly broadened. Vertical lines in Fig |2] mark 
T/v and Tg under conditions of zero stress. 

As found previously for the undoped parent com- 
pounds and for Ba(Fei_a;Co2,)2As2i'2^"— i^, pb > Pa for 
all underdoped compositions. The difference begins grad- 
ually at a temperature well above Ts, but there is no 
indication in either the resistivity or its derivatives of 
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FIG. 2: (Color online) Temperature dependence of the in-plane resistivity pa (green curves) and (red curves) of stressed 
crystals of (a) Ba(Fei_a;Nia:)2As2 and (b) Ba(Fei_i,Cui,)2As2. Data have been normalized by the value at 300 K. Values of x are 
label-led in each panel. Vertical lines mark Ts (dot-dashed line) and Tat (dashed line) determined from dp{T)/dT for unstressed 
conditions. 



an additional phase transition marking the onset of this 
behavior. The temperature at which the difference be- 
comes discernible depends on pressure^, and the effect 
appears to be associated with a large Ising nematic sus- 
ceptibility. For both series, ph rapidly develops a steep 
upturn with decreasing temperature as the dopant con- 
centration is increased from zero, similar to the behavior 
observed for Ba(Fei_3,CO:c)2As2. However, although pa 
in Ba(Fei„2:Coa;)2As2 exhibits metallic behavior up to 
X — 0.035, Pa starts increasing with decreasing tempera- 
ture in both Ba(Fei_a;Ni2:)2As2 (beginning &i x ^ 0.018) 
and in Ba(Fei_2;Cu2;)2As2 (beginning at x ^ 0.012). 
This behavior is unlikely to be associated with partial 
mixing of the 6-axis resistivity (for instance due to incom- 
plete detwinning) because pa is found to increase with 
decreasing temperature even above Ts for higher Ni and 
Cu concentrations (i.e. at a temperature for which there 
is no twin formation). 

The in-plane resistivity anisotropy, expressed as pb/ Pa, 
is shown as a function of temperature and composition 
in Figini A linear interpolation between data points 
has been used to generate the color scale images. Data 
points indicate values of Ts, Tat and Tc under condi- 
tions of zero applied stress. For comparison, data for 
Ba(Fei_a;Coa;)2As2 taken from ref. HI, are also shown. In 
all three cases, Pb/ Pa is found to vary non-monotonically 
with increasing amounts of Co, Ni or Cu. For the case of 
Co substitution, the in-plane anisotropy peaks at a value 
of nearly 2 for a composition 0.025 < a; < 0.045, close to 
the onset of the superconducting dome. Uncertainty in 
the exact composition at which Pb/ Pa is maximal reflects 
the relatively sparse data density. In comparison, for Ni 
substitution, pb/ pa peaks for 0.012 < a: < 0.022, approx- 
imately half the dopant concentration as for Co substitu- 
tion. In addition, a much weaker secondary maximum is 
found for Ba(Fei_a:Nia;)2As2 centered at a; '-^ 0.03. Mea- 
surements of multiple crystals confirmed the presence of 



this feature. For Cu substitution, the in-plane anisotropy 
peaks in the range 0.022 < a; < 0.03. 

To investigate the origin of the non-monotonic doping 
dependence of the in-plane resistivity anisotropy, which is 
observed for all three dopants, we turn now to the results 
of magnetotransport measurements, starting with the 
Hall effect (Fig|4]). For the parent compound BaFe2As2 
, it has been well established that even in modest mag- 
netic fields the transverse resistivity, pxy is non-linear22,. 
Non-linearity in p^y is expected in multi-band systems 
in which at least one FS pocket is not in the weak field 
limit. 

Addition of either Co, Ni or Cu (Figgja), (b) and (c) 
respectively) rapidly suppresses the non-linear behavior 
of Pxy The rate at which the non-linear field dependence 
of Pxy is suppressed with x is best seen by considering 
Pxy/B as a function of x for different values of i3. If pxy 
varies linearly with B, Pxy/B yields a constant value, 
which is just the Hall coefficient, Rh- Data for all three 
series are shown in FiglS]as a function of x for T = 25 K, 
where Pxy/B has been evaluated for representative fields 
5 = 0,2,5,9 and 14 T. (Data for B = were evaluated 
by considering the instantaneous slope of pxy at B — 0.) 
For the available field range (0^14 T), Pxy/B becomes 
independent of field for x > 0.045, 0.025 and 0.018 for 
Co, Ni and Cu substitution respectively. Clearly Cu sub- 
stitution is more effective at suppressing the non-linear 
Hall behavior than Ni substitution, and Ni substitution 
is more effective than Co substitution. 

It is instructive to compare the composition depen- 
dence of Pxy/B with that of the in-plane resistivity 
anisotropy, pb/pa- Black data points in Fig l5] show Pb/pa 
(referenced to the right-hand axis) also evaluated at 25 K 
as a function of x for all three substitution series. Com- 
parison of data for Co and Ni substitution reveal that 
suppression of the non-linear behavior of pxy appears to 
be correlated with the onset of the large in-plane resis- 
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FIG. 3: (Color online) In-plane resistivity anisotropy, 
pb/Pa, as a function of temperature and doping for 
(a) Ba(Fei-:,Cox)2As2, (b) Ba(Fei_^Ni^)2As2, and (c) 
Ba(Fei_a;Cux)2As2. The color scale has been obtained by a 
linear interpolation between adjacent data points. The same 
scale has been used for all three panels. Black circles, Squares, 
and triangles indicate Ts, Tn, and Tc respectively, deter- 
mined for unstressed conditions. Tc for Ba(Fei_a;Coa;)2As2 
and Ba(Fei_2:Nia;)2As2 was defined by the midpoint of the su- 
perconducting transitions, while it represents the onset tem- 
perature of the transition for Ba(Fei_a;Cu2,)2 As2. 
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FIG. 4: (Color online) Field dependence of the transverse 
resistivity p^y at a temperature T — 25K with _B || c 
for (a) Ba(Fei_xCo^)2As2, (b) Ba(Fei-xNix)2As2 and (c) 
Ba(Fei_a;Cua;)2As2. Compositions are labeled in the legends. 



tivity anisotropy. That is, Ni substitution seems to be 
almost twice as effective at both suppressing the non- 
linear Hall effect, and also at yielding a large in-plane 
anisotropy, than is Co substitution. This apparent corre- 
lation is less pronounced for the case of Cu substitution, 
for which, despite an even more rapid suppression of the 
non-linear Hall effect, the resistivity anisotropy appears 
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FIG. 5: (Color online) Left axes: composition-dependence 
of pxyjB at 25 K for (a) Ba(Fei_^Coi,)2 As2, (b) 
Ba(Fei_a;Nii,)2As2, and (C) Ba(Fei_2,Cua;)2As2. Values of 
PxyjB were evaluated at 0, 2, 5, 9 and 14 T, as described in 
the main text. Right axes (black data points): composition 
dependence of pbj po. at 25 K. 



to peak at a comparable range of compositions as found 
for Ni substitution. 

To shed more light on the correlation of the non-linear 
Hall effect and in-plane resistivity anisotropy, the trans- 
verse (B parallel to c-axis) magnetoresistance (MR, de- 
fined as IS.pl p) has also been measured on both detwinned 



and twinned samples. We first discuss the parent com- 
pound, for which representative MR data are shown in 
FigEJa) for a temperature of 25 K. As has been previ- 
ously observed, the MR of the parent compound is linear 
over a wide field range'^^. This behavior extends to very 
low fields, at which point the MR naturally reduces to a 
weak-field quadratic dependence. Interestingly, the lin- 
ear behavior does not depend strongly on the current 
direction, and the difference of the linear slope can be 
mainly accounted for by the difference of the zero field 
resistivity, i.e. Ap/p scales approximately with B j p. As 
was first shown by Abrikosov, a linear band dispersion 
can lead to a linear MR in the quantum limit2^"=2I. In 
the case of BaFe2As2, the linear MR could be naturally 
explained by the presence of Dirac pockets in the AFM 
reconstructed state due to the symmetry protected band 
crossing^. Since the Dirac pockets have a small volume 
and a long mean free path-'--, the high field limit can be 
reached with moderate fields, which could account for the 
non-linear Hall effect observed for the parent compound 
described above. 

The cross-over from the weak-field dependence to 
the high-field linear dependence can be best seen by con- 
sidering the field derivative of the MR, (c?[Ap/p]/c?B), 
which is plotted in FigE] for the parent compound. At 
low fields, Ap/p = resulting in a linear field de- 

pendence for d\^pl p\l dB as B approaches zero. How- 
ever, above a characteristic field B*, d\/S.p/ p]/ dB starts 
to deviate from this weak field behavior, and appears 
to saturate to a much reduced slope. This indicates 
that at high field the MR is dominated by a linear field 
dependence, but there is also a small quadratic term 
{Aplp^ArB + 0{B^)). 

Substitution of Co or Ni rapidly suppresses the linear 
MR observed for the parent compound. Representative 
data are shown in Figl6l^b) for a detwinned single crystal 
of Ba(Fei_2:Coa;)2As2 with x — 0.035 at 25 K. As can 
be seen, a linear MR is still observed, but the weak-field 
quadratic behavior extends to a higher field values. As 
for the parent compound, the anisotropy in the linear 
slope can be mainly accounted for by the anisotropy in 
Pa and pfc. 

A comprehensive doping-dependence was obtained 
for twinned samples of Ba(Fei_j;Coa;)2As2 and 
Ba(Fei_a;Nia;)2As2. Representative data are shown 
in Figs. [S] (d) and (e) respectively, illustrating the 
rapid suppression of the MR with substitution. As 
described for the parent compound, we can extract 
the characteristic field B* by considering the field 
derivative of the Ap/p (Fig. [T]). Fitting the high 
field [B > B*) MR by a second order polynomial, we 
obtain the coefficient of the linear field dependence 
coefficient Ai, the doping dependence of which is shown 
in Fig. [71 As can be seen, the characteristic field scale 
B* increases rapidly as a function of doping, whereas 
the linear coefficient decreases and almost vanishes at 
X = 0.051 for Ba(Fei_^C0:r)2As2 and x = 0.027 for 
Ba(Fei_a;Co2;)2As2. Apparently, the rate at which the 
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FIG. 6: (Color online) Transverse magnetoresistance(Mii) 
(Ap/p) of Ba(Fei_^Co^)2As2 for (a) a; = (b) 3.5% at 
T = 2'b K. The MR has been measured on detwinned sam- 
ples with current along a/6 axis (Green/Red curves), and also 
on twinned sample(Black curve). (c,d) The MR of twinned 
samples of Ba(Fei_3;Co2,)2As2 and Ba(Fei_a;Nia;)2As2 respec- 
tively at 25 K. *Data for a; = has been scaled down by a 
factor of two for clarity. 



linear MR is suppressed is twice as rapid for Ni substi- 
tution as for Co substitution. This behavior is clearly 
correlated with the suppression of non-linearity in Hall 
effect, which occurs over a similar range of compositions 
as shown in Fig. [SJ providing additional evidence that 
the non-linear Hall coefficient in the parent compound is 
due to the presence of high-mobility Dirac pocket(s). 



IV. DISCUSSION 

All three substitutions investigated in this study ex- 
hibit a non-monotonic variation of the in-plane resistiv- 
ity anisotropy as the dopant concentration is progres- 
sively increased (Fig. Without further information 
it is not clear whether this effect is related to changes 
in the anisotropy of the scattering rate, or to changes in 
the electronic structure. However, as described above, 
consideration of the magnetotransport properties is sug- 
gestive of an important role for the Dirac pocket of the 
reconstructed Fermi surface. The progressive suppres- 
sion of the linear MR with chemical substitution, and 
the associated suppression of the non-linear Hall coef- 
ficient, point to a scenario in which the contribution to 
the conductivity from the FS pockets associated with the 
protected band crossing (the Dirac pockets) rapidly di- 
minishes with increasing dopant concentrations. For a 
multi-band system, the conductivity tensor is the sum 
of the contribution from each Fermi surface. If one par- 
ticular Fermi surface pocket dominates the conductivity 
tensor, then the transport anisotropy will also be deter- 
mined by the anisotropy of that particular Fermi surface. 
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FIG. 7: (Color online) (a) The field derivative of 
MR of twinned samples of Ba(Fei_a;Coa;)2As2 and (b) 
Ba(Fei_a;Ni2,)2As2 *Data for x = Q has been scaled down 
by a factor of two for clarity, (c) The field derivative of MR 
(dMR/dB) measured on a twinned sample of the parent com- 
pound at r = 25K. A critical field scale B* clearly divides 
the MR behavior into two regimes: below B* the MR shows 
a weak field quadratic behavior and above B* MR shows 
a high field linear behavior, (d) The doping evolution of 
the field scale B* (blue circles) and the high field MR lin- 
ear coefficient Ai (red squares), for Ba(Fei_i:Coi:)2As2 and 
Ba(Fei_a;Nia;)2As2 (solid and open symbols respectively). 
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As observed previously by photoemission measurements, 
the Dirac pockets have an almost isotropic in-plane Fermi 
velocity^. If the mobility of this pocket is such that it 
dominates the transport, it would severely diminish the 
anisotropy associated with any other Fermi surfaces, just 
as we observe for the parent compound. The contribution 
from these Dirac pockets is progressively weakened by the 
transition metal substitution, which is manifested in the 
magnetoresistance and Hall effect. The subsequent emer- 
gence of a large in-plane resistivity anisotropy clearly in- 
dicates that the remaining low-mobility FS pockets are 
highly anisotropic. This change is also manifested in the 
normalized value of the residual resistivity at low tem- 
perature, which shows an abrupt increase at the same 
composition (Fig. [9]). 

The mechanism that suppresses the contribution from 
the Dirac pockets is unclear. The band crossing is pro- 
tected by crystal inversion symmetry, but introducing 
impurities into FeAs planes locally breaks this symme- 
try. This effect would not only open a gap at the Dirac 
point, but would also increase the scattering rate for the 
Dirac electrons due to mixing of the orbital wavefunc- 
tions. This is consistent with the reduced mobility de- 
rived from the MR analysis, and also consistent with the 
suppressed non-linearity in the Hall effect. On the other 
hand, Co substitution is argued to effectively electron- 
dope the system, which shifts the chemical potential. 
If there is a gap at the Dirac point, this chemical po- 
tential shift could possibly lead to a Lifshitz transition 
in which the Dirac pocket vanishes. This possibility has 
been extensively discussed in the recent papers by Liu et 
al J^iii A recent Nernst effect measurement also shows a 
suppression of Dirac transport by Co substitution in the 
Eu(Fei_;rCO;r)2As2 systeirt^, suggesting the effect is not 
restricted to BaFe2As2 system. 

To make a more quantitative understanding of how 
the Dirac pockets are being suppressed as a function of 
doping, one would ideally like to obtain the transport 
parameters for each band. According to Local Density 
Approximation (LDA) calculation, the BaFe2As2 parent 
compound has four closed Fermi surface pockets in the 
AF reconstructed states, therefore one would need eight 
independent parameters (mobility and carrier density for 
each pocket) to characterize the transport properties. If 
we focus only on the data for fields close to zero, we 
can obtain the coefficient of the quadratic term, A2, 
by inverting the sum of the conductivity tensors of the 
four Fermi surfaces. We denote the conductivity and 
mobility of each Fermi surface by aij and iiij , where the 
index i — e,h stands for electrons or holes, and the index 
j — D,P represents the Dirac bands and the parabolic 
bands.. We assume no intrinsic magnetoresistance for 
each individual Fermi surface, because in the condition of 
isotropic scattering rate and an ellipsoidal FS the leading 
quadratic term in Zener- Jones expansion is zero. The 
parameter A2 is given by 



A. 



((Te + (JhY 



-A2,e + - 



-A2,h (1) 



A2.^ = 



0'i,pO'i,D(Mi,_D — ^^i,pY 



Mi = ; ^^i,D 



-^^i,P 



O'i^D + <^i.P- 



(2) 



(3) 



(4) 



Here, a^, ah, fJ-e, fJ-h are the effective electron and hole 
conductivity and mobility in zero field. If we assume 
the the Dirac bands are dominating the transport, i.e. 

Cre,D,Crh,D ^ ^e,PiCh,P and He,D,fJ-h,D ^ fJ'e,P , fJ-h,D , 

the above expression can be greatly simplified because 
^ <yi,D and ~ /^i,_D. The second and third term in 
equation 1 can be assumed to be much smaller than the 
first term: 



-Ao 



(cTe + (Th){cre,D + CTe.p)^ 



0'e,PO'e,D(A*e,P — Mc.d)^ 



CTe.PMe 



0'e,p(o'e + <^h)^^ 



2 a^ahifJ-e + IJ-hY 



((Te + Crft)2 



(5) 
(6) 
(7) 



and equation 1 reduces to only the first term, which only 
depends on the effective mobility and conductivity of 
electrons and holes. The simplicity of this expression al- 
low us to make a physical interpretation of the coefficient 
A2 in terms of an effective mobility. 



■ (^e + Ai/i) = [i-MR < -(He+ fJ-hj 



CTh 



.(8) 



The square root of the quadratic field coefficient A2, 
which we denote as ^mr, gives the lower bound of the 
average mobility of electrons and holes (^aue)- The equal- 
ity between ^ave and ^mr holds only when the electron 
and hole conductivities are equal. However, even in the 
case of a strongly asymmetric conduction scenario, the 
Mmp still gives a good estimate of ^ave- For example if 
fe/c/i = 10 then ^mr = O.G/Ua^e- AH of the above rea- 
soning makes sense only when the assumption that the 
Dirac bands dominate the transport holds. To be consis- 
tent, we extract the real number of the parent compounds 
to see if this is really the case. The value of fi mr of the 
parent compounds at 25 K is about 1130 cm^/Vs, which 
is comparable to the mobility of the Dirac pockets ex- 
tracted from quantum oscillations for crystals prepared 
under similar conditions (^ 1000 cm^/Vs)®. This is in- 
deed consistent with the assumption of /ig ^ Me. pi- We 
can also obtain the effective carrier density umr. by us- 
ing the relation a = ne/i. If (Te ~ (Te,£), then umr. should 
be similar to the carrier density of Dirac carriers rather 
than the total carrier density. In the parent compound 
the observed value of (Xmr = 1130 cm^/Vs corresponds 
to a very low effective carrier concentration, umr = 0.003 
electron per Fe. This value is much lower than the total 
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number of carriers that one would obtain from the LDA 
calculations and ARPES measurements, but is compara- 
ble to the size of Dirac pockets as observed from the quan- 
tum oscillation experiments^. Therefore the assumption 

Ce O'e.D slsO holds. 

The effective mobility extracted from MR is much 
higher than the mobility directly obtained from the Hall 
coefBcient fiHaii = l-Ro"! = 376cm^/ys. This is because 
the contribution of electron and hole mobility in the Hall 
coefficient cancel each other: 



R 



at Re 



fJ-Hall 



- ^hfJ-h 



(9) 
(10) 



In fact, the much smaller value of Hall mobility than 
the MR mobility already implies that the contribution to 
the conductivity of electrons and holes are of the same 
scale. To understand this one can consider the opposite 
case, one where electrons dominate the conduction /Xe 3> 
fih and CTe > o-^: 



(TeA*e + CT/iM/i 










C^e + at 



f^e (11) 



-A*e (12) 



< 1 fJ.AIR < ^Hall (13) 

Therefore our measured value suggests that both elec- 
trons and holes play an important role in the transport 
in the reconstructed state in the parent compound, and 
this also gives us the confidence that ^j.mr is a good es- 
timate of fJ-ave, since their difference is smaller when Ue 
gets closer to ah- 

As we argued above: the high (Imr and small hmr. re- 
flects the fact that transport in the parent compound is 
dominated by a small number of high mobility carriers, 
i.e. the carriers from the Dirac pockets. However, this 
is no longer the case as we increase the Ni or Co doping 
concentration. The extracted value of (J-mr and umr as 
a function of doping is plotted in Fig. [51 By increasing 
doping concentration fiMR decreases rapidly and umr 
increases rapidly. Apparently our previous assumption 
that helped us simplify the MR expression is no longer 
valid, and it is difficult to make a simple physical inter- 
pretation of the extracted fiMR or umr- Nevertheless, 
the observed doping evolution is highly suggestive of a 
shift of the dominant role in transport to the high car- 
rier density and low mobility carriers. Again the effect 
of Ni doping on suppressing the Dirac carriers is twice 
as fast as Co. Our data reveal that as the conductiv- 
ity from the Dirac pockets is progressively suppressed, 
a large in-plane resistivity anisotropy emerges. The di- 
rect implication is that the other pockets of reconstructed 
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FIG. 8: (Color online) The doping evolution of the effective 
MR mobility (imr and effective AIR carrier density umr 
of Ba(Fei_2:Cox)2As2 and Ba(Fei_2:Nia;)2As2 extracted from 
weak field MR, as described in the main text. 



FS are highly anisotropic, which is borne out by recent 
quantum oscillation measurements of BaFe2As2 3 

Evidence for the suppression of the contribution of a 
high mobility pocket of reconstructed FS can also be 
found in the doping-dependence of the resistivity at low 
temperature. A direct comparison of the magnitude of 
the in-plane resistivity normalized by its room tempera- 
ture value^^^ and the in-plane resistivity anisotropy at 
T = 25 K is plotted if Fig. [SI As can be seen, the 
two quantities follow each other closely in the case of Co 
and Ni doping, but not for Cu substitution. Based on 
the previous analysis, the doping evolution of the resis- 
tivity can be readily understood. Initial suppression of 
the contribution to the transport arising from the Dirac 
pocket leads to an initial rise of the normalized resistiv- 
ity. With progressive doping the magnetic order is fur- 
ther suppressed, releasing carriers and hence leading to 
a decrease in the resistivity for higher dopant concentra- 
tions. Note that the progressive doping also suppresses 
the structural transition that breaks the rotational sym- 
metry, therefore after reaching a maximum value, the 
in-plane anisotropy also decreases. 

We now comment briefly on the trends revealed by 
comparison of Co, Ni and Cu substitution. First princi- 
ples density functional calculations indicate that the ad- 
ditional charge associated with Co and Ni substitution in 
BaFe2As2 resides within the muffin-tin potential*^, but 
an associated rigid-band shift leads to an effect that is ul- 
timately equivalent to electron doping^^. In calculations 
for transition metal impurities in LaFeAsO, Ni is found 
to be approximately twice as effective as Co in terms of 
the rigid band shift^, motivating a direct comparison 
of results for Ba(Fei_a;Coa:)2As2 and Ba(Fei_j;Ni2;)2As2. 
In contrast, similar calculations for Zn and Cu impuri- 
ties appear to strongly affect the electronic structure^i^. 
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and (c)Ba(Fei_i,Cui,)2As2. 



possibly accounting for significant differences between 
the phase diagrams of Ba(Fei_a;Cua;)2As2 and those of 
Co and Ni substituted BaFe2As2, including the much 
weaker superconductivity in Cu-substituted BaFe2As2^^. 
As anticipated, Ni substitution suppresses the contribu- 
tion of Dirac carriers to the transport. Furthermore, this 
suppression occurs twice as rapidly with x for Ni as for 
Co substitution. Similarly, the onset of the large in-plane 
resistivity anisotropy is found to occur for approximately 
half of the values of x than is the case for Co substitu- 
tion (Fig|3]). The origin of the secondary maximum in 
the in-plane anisotropy observed for Ba(Fei_a;Nia;)2As2 
is unclear, but perhaps reflects subtle changes in the re- 
constructed FS with chemical substitution. 

Inspection of FigO reveals that Cu substitution sup- 
presses the non-linear Hall coefficient even more rapidly 
than Ni substitution. This effect is consistent with the 
deeper impurity potential of Cu relative to Ni, acting 
to increase the elastic scattering rate more rapidly per 



impurity added. However, the onset of the large in- 
plane resistivity anisotropy occurs for somewhat larger 
compositions, peaking at x ^ 0.024. Given that the 
effects of band filling and scattering are believed to be 
somewhat different for Cu substitution relative to Co 
or Ni substitution, it is perhaps not surprising to find 
that the transport properties evolve in a slightly differ- 
ent manner. Nevertheless, it is interesting to note that 
despite this, a region of the phase diagram still exists over 
which a large in-plane resistivity anisotropy is observed, 
presumably because high mobility isotropic FS pockets 
have been suppressed. Perhaps significantly, the compo- 
sitional range over which this anisotropy is observed is 
narrower for Cu substitution than it is for Ni substitu- 
tion, which in turn is somewhat narrower than for Co 
substitution. 

Finally, an intriguing correspondence can be made be- 
tween the iron pnictides and underdoped cuprates42i Re- 
cent Nernst measurements reveal a large in-plane elec- 
tronic anisotropy onsets at the pseudogap temperature^. 
Further Hall coefficient and quantum oscillation mea- 
surements suggested at even lower temperatures broken 
translational symmetry causes a reconstruction the Fermi 
surface and a high mobility isotropic electron pocket 
emerges^. This high mobility electron pocket dominates 
the low temperature transport. As a result, not only 
does the Hall effect change sign from positive to neg- 
ative, but the large electronic anisotropy which onsets 
at pseudogap temperature is also reduced. It has also 
been suggested that at a critical doping the system could 
undergo a Lifshitz transition at which the high mobility 
electron pocket disappears, which is also accompanied by 
an enhancement of in-plane resistivity anisotropy^. 



V. CONCLUSIONS 

Co, Ni and Cu substituted BaFe2As2 are all found 
to exhibit a large in-plane resistivity anisotropy over 
a certain compositional range on the underdoped side 
of the phase diagram. The non-monotonic variation in 
the resistivity anisotropy as the dopant concentration is 
increased is especially striking given the uniform sup- 
pression of the lattice orthorhombicity. The non-linear 
Hall coefficient and linear MR which are observed for 
the parent compound^i^, and which are likely associ- 
ated with the Dirac pockets of the reconstructed FS, are 
suppressed with increasing dopant concentrations. In- 
triguingly, for both Co and Ni substitution, for which a 
direct comparison is motivated based both on their re- 
spective phase diagrams and also on first principles LDA 
calculations—, the large in-plane resistivity anisotropy is 
found to emerge over the same range of compositions 
at which the non-linear Hall and linear Mi? are pro- 
gressively suppressed. Consideration of this evidence 
suggests that the isotropic, high- mobility Dirac pock- 
ets revealed by dHvAs^'^'^ ARPES^ and magnetotrans- 
port measurements^, might effectively mask the intrinsic 
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in-plane transport anisotropy associated with the other 
pockets of reconstructed FS. Within such a scenario, only 
when the contribution to the conductivity from the Dirac 
pockets is suppressed can the underlying anisotropy be 
revealed in the transport, perhaps accounting for the 
non-monotonic doping dependence. 

Finally, we note that it remains to be seen to what 
extent the presence of a large resistivity anisotropy is a 
generic feature of the phase diagram of Fe-pnictide su- 
perconductors. Recent measurements of K-substituted 
BaFe2As2 indicate that the hole doped analog may not 
exhibit a large in-plane anisotropj*2i,. This might reflect 
differences in the effect of electron vs hole doping on the 
reconstructed FS, or perhaps differences in the elastic 
scattering rate, since chemical substitution away from 



the Fe plane will presumably have a weaker effect. 
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